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ABSTRACT 

This  progress  report  covers  the  first  15  months  of  a  continuing  study 
on  construction  of  a  comprehensive  mathematical  model  of  composite  solid 
propellant  combustion,  with  emphasis  on  work  during  the  past  12  months 
since  reporting  on  tiie  early  efforts  at  the  meeting  last  year. 

The  liquid  layer  believed  to  exist  on  burning  ammonium  perchlorate 
has  been  treated  in  a  manner  similar  to  that  described  last  year.  Depth 
of  penetration  of  a  thermal  wave  is  calculated  and  this  value  is  multiplied 
by  the  surface  area  of  oxidizer  on  the  burning  propellant  surface.  The 
surface  is  calculated  in  the  computer  program  for  use  in  other  equations. 
Totni  volume  of  the  liquid  layer  derived  in  this  manner  is  then  used  to 
calculate  the  thermal  contribution  to  surfnee  temperature. 


Treatment  of  aluminum  has  been  added,  and  the  treatment  of  multimodal 
oxidizer  has  been  improved.  Most  internal  parameters  have  been  thoroughly 
checked  to  verify  that  they  have  realistic  vnlues.  Considerable  time  was 
spent  on  the  Burke-Schumann  flame  height  calculations.  It  has  been  found 
that  at  least  seven  terms  of  the  series  solution  must  bo  considered.  The 
single  first-term  calculation:?  that  have  been  used  in  the  pa s t  do  not  give 
a  good  assessment  of  flame  height:  the  flame  height  is  expressed  implicitly 
in  the  final  equation,  and  it  was  necessary  to  solve  it  by  a  Newton-Rnphson 
iteration.  Improvement  in  a  number  of  input  parameters  is  described.  Gas 
phase  conductivities  were  calculated  using  Cha pmnn-Ensc og  theories  with 
the  llirschfelder-Dicken  corrector.  The  same  theories  were  used  to  calculate 
diffusion  coefficients.  Stoichiometric  ratios  were  calculated  using  more 
realistic  pyrolysis  species.  A  thermogravimetrie  analysis  was  made  of  a 
polybutadiene  propellant  binder  and  a  statistical  method  was  used  to  derive 
kinetic  parameters  from  the  thermograms.  The  method  used  was  the  Zsake 
modification  of  Doyle's  method,  which  is  probnbiy  the  most  sophisticated 
way  of  treating  TGA  data.  Some  other  improvements  in  programming  and  cor¬ 
relation  with  actual  propellant  motor  firings  are  shown  and  future  plans 
are  given. 


INTRODUCTI ON 


Work  continued  during  the  past  year  on  the  comprehensive  modeling  of 
composite  solid  propellants.*  The  model  being  developed  is  based  on  the 
Becksteod,  Derr,  I’rice  multi  flame  rodel  (i,2)  and  utilizes  a  critical  par¬ 
ticle  size  theory  (“)  to  provide  n  moons  of  calculating  condensed  phase 
heat  release.  A  number  of  improvements  have  been  made  to  the  model  and  the 
computer  program  that  was  described  at  the  meeting  last  year  (;i);  one  of 
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the  major  improvements  being  the  Addition  of  a  treatment  of  the  liquid 
layer  an  burning  ammonium  perchlorate  (AP).  Tiie.se  improvements,  .mi  exam¬ 
ple  experiment  and  future  plans  are  described  in  the  following  paragraphs. 

MATIltMAflCAL  MODEL  \ND  IMPHW'KMKCTS 

The  most  recent  surface  energy  balance  is  shown  on  Fig.  1.  The  new 
terms  (not  shown  last  year)  are  the  energy  to  melt  the  aluminum  and  the 
energy  from  the  liquid  layer.  In  this  equation  mox  is  the  mass  flux  of 
interstitial  propellant;  or,  more  specifically,  the  hinder  plus  the  sub- 
critical  AP.  The  term  for  condensed  phase  heat  release  based  on  the 
critical  particle  size  theory  is  included.  The  critical  particle  size 
theory,  stated  simply,  is  that  there  exists  for  a  given  propellant  and 
pressure  a  particle  size  where  all  smaller  particles  will  be  consumed  in 
condensed  phns  ^  reaction.  The  terms  for  aluminum  fusion  and  liquid  layer 
on  AP  are  also  present.  The  treatment  of  aluminum  used  is  somewhat  dif¬ 
ferent  from  that  used  in  the  BDP  model;  with  the  methodology  explained 
later. 


Fig.  1.  Critical  Particle  Size  Model  Energy  Balance 


Figure  2  shows  first  the  current  equation  for  surface  temperature. 

This  was  derived  from  the  energy  equation  using  mass  balance  equations  in 
the  same  manner  as  orininnllv  used  by  Becks  tend  (2).  The  next  three  equa¬ 
tions  are  for  the  heat  release  in  the  primary,  final,  and  AP  flames. 

Flame  Heights.  The  Burke-Schumnnn  equation  modified  for  short  flames  is 
shown  in  Fig.  3-  In  the  past  a  single  term  of  this  expression  has  been 
used  but  we  have  found  that  at  least  seven  terms  should  be  used.  The  equa¬ 
tion  is  really  not  ns  complicated  ns  it  first  appear*:  it  is  nothing  more 
than  is  shown  in  the  second  equation.  By  adding  the  extra  terms  the  equa¬ 
tion  has  been  made  implicit  in  flame  height  but  it  is  solved  easily  by  a 
Newton-Ihiphson  iteration.  Convergence  is  obtained  rapidly  when  the  initial 
guess  for  flame  height  is  taken  ns  one-tenth  the  value  obtained  from  a 
single-term  calculation.  This  table  is  the  numerical  value  of  the  first 
five  terms;  and  each  term  is  about  hair  the  value  of  the  preceding  term. 

An  actual  check  of  flame  heights  when  different  numbers  of  terms  are  used 
shows  that  at  least  seven  terms  must  he  used  to  get  the  flame  within  1 
micron. 
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Fig.  2.  Current  Equations  for  Surface  Temperature 
Used  in  Computer  Program 
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Fig.  3>  Burke-Schuma tin  Equa ti on  Modified 
for  Short  Flames 
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A  number  of  other  improvements  have  been  made  to  the  computer  pro¬ 
gram  for  this  model.  lXiring  the  past  year  nil  parameters  have  been 
thoroughly  checked,  particularly  those  used  internally  and  normally  not 
listed.  We  found  that  some  of  the  volume  and  weight  fractions  would  be¬ 
come  unrealistic;  e.g.,  greater  than  one.  When  this  occurred  it  caused 
discontinuities  in  the  burn  rate  curve.  This  was  remedied  by  setting 
boundaries  on  these  parameters. 

Treatment  of  the  oxidizer  particle  size  has  been  improved.  The  spline 
fit  of  the  particle  size  distribution  curve  is  still  used.  A  normal  hi-  or 
tri-modal  oxidizer  blond  will  usually  be  split  into  10  modes  by  this  tech¬ 
nique.  Calculation  of  the  h/l)  ratios  is  handled  in  a  different  manner 
than  that  described  by  Cohen,  Derr,  and  Price  last  year  (Fig.  4).  Each 
mode  of  supercritical  oxidizer  is  split  into  50  parts;  the  ignition  delay 
term  in  the  equation  for  h/D  is  calculated  for  each  part,  summed  and 
averaged.  Each  ignition  term  is  then  multiplied  by  the  weight  fraction 
of  the  mode  and  summed  with  similar  terms  produced  from  each  mode.  Then 
h/D  and  Sox  are  calculated. 


Fig.  4.  Calculations  for  h/D  Ratios 

Oxidi zer/fuel  stoichiometric  ratios  are  different  in  the  primary  flame 
and  the  final  flame.  At  temperatures  such  ns  are  found  at  the  binder  sur¬ 
face  it  is  known  that  polybutndiene  depolymorizes  (5)  with  the  major  specie 
produced  being  butadiene.  This  is  considered  the  major  diffusing  specie 
in  the  primary  flame,  giving  a  stoichiometric  ratio  of  3. 7-  The  Chapmnn- 
Encog  equation  in  the  center  of  Fig.  4  was  then  used  to  calculate  the  dif¬ 
fusion  coefficient  for  the  primary  flame.  In  the  final  flame  it  was 
assumed  that  hydrocarbon  fragments  were  the  major  diffusing  species;  a 
stoichiometric  ratio  of  8.0  was  used,  along  with  a  larger  diffusion  co¬ 
efficient. 

Liquid  Layer.  The  AP  liquid  layer  was  not  treated  ns  a  melt  but  ns  an 
exothermic  layer.  We  have  always  assumed  this  layer  to  be  a  solution  of 
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AP  and  perchloric  acid.  A  thermal  wave  penetration  is  calculated  m  the 
same  manner  as  for  the  critical  particle  size.  This  volume  is  converted 
to  weight  fraction  of  A1’  in  the  propellant  that  is  available  for  exothermic 
surface  reaction  (Kef  equation  at  bottom  of  l’ig.  4). 


EXPERIMENTAL  RESULTS 


The  combustion  model  with  the  improvements  noted  above  was  applied  to 
several  smokeless  propellants  for  which  experimental  data  were  available. 
Theoretical  calculations  for  three  of  the  propellants  (based  on  carboxy- 
terrainated  polybutadiene  binders)  were  made  at  two  initial  propellant 
tempera tores  and  temperature  sensitivities  of  burning  rates  were  calculated. 
As  shown  in  Table  1,  agreement  of  77  F  burn  rates  is  good.  Slopes  of  the 
burn  rate  curves  are  identical.  The  <7^  values  were  calculated  using  the 
one  over  one  minus  "n"  relationship.  The  theoretical  Cq(  values  are  all  low 
and  should  be  two  to  three  times  higher;  and  we  will  look  into  this  further. 


Table  I.  Temperature  Sensitivity 
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Figure  5  shows  a 
burning  rate  vs  pres¬ 
sure  plot;  the  upper 
line  represents  a  least 
squares  fit  to  rocket 
motor  firings;  the  lower 
line,  the  theoretical. 

This  is  a  two  cycle -by- 
two  cycle  log-log  plot. 

The  scales  are  the  same 
as  those  widely  used  in 
industry  to  plot  rocket 
motor  rate;  therefore, 
the  da  ta  are  presented 
in  a  realistic,  prac¬ 
tical  manner.  As  noted 
last  year  we  did  not  use 
strand  data  for  corre¬ 
lative  purposes  since 
strand  rates  are  not 
necessarily  predictive 
of  motor  rates.  Pro¬ 
pellant  used  in  the  motor  from  which  the  plotted  data  were  taken  has  a 
Butarez  CTL-type  binder  and  an  8bjr  solids  loading.  The  trimodal  oxidizer 
is  comprised  of  commercial  200-  and  80-micron  fractions  and  Ma jac-ground 
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Propellant  with  Butarez 
86^  Solids  and  Trimodal 


CTI.  Type 
Oxidir  r 


wi  th 


10-micron  Al'  at  weight  percent  levels  35/26/25.  The  theoretical  curve  is 
identical  in  slope  up  to  1000  psi. 


Comparable  data  0.8 
with  a  propellant  con¬ 
taining  a  trimodal  0 6 

mixture  of  400- ,  200- , 
and  80-tuicron  AP  at  05 

weight  percent  levels  u 
51/20/15,  respectively,  S*  '  { 
are  shown  in  Fig.  6.  * 

The  break  at  1000  psi  w-03 
is  well  verified  by  g 
motor  firings.  No  g 
curve  fitting  was  used  SO. 2 
to  obtain  these  re¬ 
sults;  i.e.,  they  were 
obtained  on  the  first 
computer  run.  Fig.  7 
shows  data  from  a  0  , 

hydroxy- terminated  ? 

polybutadiene  (HTPB)- 
based  propellant  wi th 
a  large  amount  of  fine 
oxidizer;  i.e.,  a  tri¬ 
modal  blend  of  200-, 

80- ,  and  10-micron  at 
weight  percent  levels 
of  Jl/l5/%0,  respec¬ 
tively.  Again,  we  are 
comparing  motor  data; 
and,  although  the  burn 
rate  is  somewhat  lower, 
the  two  lines  have 
identical  slopes.  Data 
from  a  propellant  that 
has  a  tetramodal  oxi¬ 
dizer,  i.e..  400- .  200- , 
80-,  and  10-micron  at 
weight  percent  levels 
of  56/13/12/5,  respec¬ 
tively,  are  shown  in 
Fig.  8.  We  would  like 
to  emphasize  again  that 
no  manipulation  of 
parameters  was  required; 
the  curve  represents  the 
firslt  run  made  for  the 
oxidizer  blend  with  the 
current  computer  program. 
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Fig.  6.  Propellant  with  Butarez  CTF  Type 
Binder  with  86$  Solids  and 
Trimodal  Oxidizer 
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Fig.  7.  Propellant  with  HTPB-Type  Binder 
with  86$  Solids  and  Trimodal 
Oxidizer 
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Fig.  8.  Propellant  with  HTPB-Type  Binder 
with  86$  Solids  and  Trimodal 
Oxidizer 

FUTURE  PLANS 

In  the  ensuing  months,  wo  plan  to  further  study  the  possibilities  of 
a  two  surface  temperature  model,  add  a  treatment  for  the  catalyst,  and 
continue  studying  the  effects  of  particle  size.  We  will  also  establish 
possible  maximum  and  minimum  burn  rales  and  continue  to  make  improvements 
in  the  computer  program. 
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